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Abstract 

We reinvestigate the leading QCD corrections to the radiative de- 
cay b ^ for supersymmetric extensions of the Standard Model. 
Although the major contributions to the corrections originate from 
the running of the effective Lagrangian from the W scale down to the 
b scale, additional corrections are expected from large mass splittings 
between the particles running in the loops, as well as from integrating 
out heavy particles at scales different from the W mass. The calcula- 
tion is performed in the framework of effective field theories. 
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1 Introduction 



Among the rare decays of B mesons, the recently observed radiative weak 
decays B — ^ Xs7 [IH, where Xg is a hadronic state with total strangeness 
S" = — 1, have received much attention. As a loop-induced FCNC process it 
is in particular sensitive to contributions from physics beyond the Standard 
Model (SM). 

Since the h quark mass is much larger than the QCD scale A, one assumes 
that the inclusive decay rate is well described by the spectator model, where 
the h quark undergoes a radiative decay. The transition amplitude is given 
by the matrix element an effective magnetic moment operator. To lowest 
order 0, the coefficient of this operator is obtained by integrating out all 
heavy particles [t quark, W boson, . . . ), leaving one with an effective field 
theory describing the transition 6 — > 57 at the parton level at the weak scale. 
The QCD corrections to this coefficientB have been calculated to leading 
logarithmic accuracy in 0, ^ H, 0, |], ^ and are known to enhance the rate 
within the SM by a factor of 2-4, depending on the masses of the h and of 
the t quark. This enhancement is, however, subject to large uncertainties 
due to the poor knowledge of some input parameters like the strong coupling 
constant as and due to the residual renormalization scale dependence (for a 



recent discussion in the context of 6 ^ S7 see e.g. [^0|, |TT|, 0), which we will 
however not address in this work. 

On the other hand, if the particles in the loop have vastly different masses, 
one expects sizeable corrections to the Wilson coefficients already at the weak 
scale. These contributions, which are usually considered as a next-to-leading 
order effect, have been discussed in ref. for the Standard Model where, 
in the case of a top quark much heavier than the W, they were found to give 
an additional enhancement of the order of 20%. 

In the case of the Minimal Supersymmetric Standard Model (MSSM) 
the situation is even more complicated. First of all, due to the richer parti- 
cle content, there are more diagrams contributing to the magnetic moment 
operators, and due to the larger number of free parameters there are many 
potential additional sources of flavor changing neutral currents |15|]. How- 



ever, if one assumes further that the MSSM is a low energy effective theory 



^For a recent review and earlier work see e.g. ref. 



1 



from minimal supergravity [ITB] with radiative breaking of the electroweak 



symmetry, it is known |17, 11, a|, 19, EO, 21, 22, El, 24, 25 1 that, besides 



the SM contribution mediated by the W, there are only significant contri- 
butions from the charged Higgs (H^) and from the chargino (xi^2) exchange 
diagrams, while already with the present experimental lower limits on the 
supersymmetry (SUSY) spectrum the gluino contribution is small and the 
neutralino contribution is always negligible.i The W and the contribu- 
tions always have the same sign in the MSSM, but the chargino contribution 
can have either sign and may e.g. cancel the contribution or (for small 
chargino masses and large tan/3) even dominate the amplitude. 

Since SUSY has to be broken, the mass splitting between the various particles 
running in the loop can be very large, leading to additional important QCD 
corrections. We will advocate that, in a parameter space analysis in the 
MSSM, one should not simply add up the contributions of all diagrams at 
the W scale and use the renormalization group evolution to run down this 
sum to the b scale, but rather consider the individual contributions separately. 
This is especially important for the chargino contribution, since the lightest 
chargino can be significantly lighter than the W. 

For the reasons mentioned above, we will ignore the contributions from dia- 
grams with gluinos and neutralinos in the present work. They may easily be 
included; corrections of the type considered in this work will however always 
be numerically unimportant. 

Our strategy will be similar to the work by Cho and Grinstein |]TB[. Starting 



from the full theory at sufficiently high scales, we will construct a series of 
effective theories that is well suited for the description of the low-energy 
physics of interest. We shall give all ingredients that are necessary to obtain 
the leading QCD corrections to the 6 — > S7 inclusive rate and discuss some 
simple estimates for an MSSM scenario with the assumptions mentioned 
above. A full parameter space analysis, which depends on the details of the 
implementation of the soft SUSY-breaking, is however beyond the scope of 
the present paper and will be discussed elsewhere. 



^This assumption is supported by present experimental data on B — » X^J as well as 
the lack of evidence for large contributions beyond the SM to other FCNC processes, e.g. 
j^Oj^o niixing and rare K decays. For a discussion of a scenario with a light gluino in the 
mass range 2-5 GeV see e.g. j26| (and references quoted therein). 
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This paper is organized as follows: In section 2 we briefly review the elements 
of effective field theories needed for the present work. Section 3 explains in 
detail the calculation for a type-II two-Higgs doublet model, which is con- 
tained in the MSSM, while the contributions of SUSY particles are discussed 
in section 4. We shall present our results in section 5 and finally conclude. 



2 Effective field theory and b ^ 

The basic idea of effective field theories is by now well established, and many 



excellent reviews have appeared in the literature p7| , |28| , |29| , [30| , pl|] . Start- 
ing from some underlying full theory, one integrates out the heavy degrees 
of freedom, thereby producing a tower of non-renormalizable interactions 
(with couplings proportional to inverse powers of the heavy particle mass) 
which contain the virtual heavy particle effects. One then runs the result- 
ing effective field theory down to the appropriate scale of interest using the 
renormalization group. If additional heavy particle thresholds are crossed 
during the renormalization group running, then these particles will also be 
integrated out. The major advantages of using an effective theory for the 
calculation of low-energy observables are convenience, since calculations are 
usually simpler than in the full theory, and the gain of insight. 

A nontrivial feature of the effective field theory framework is the auto- 
matic summation of large logarithms that originate from perturbatively cal- 
culable short- distance physics by the renormalization group. As explained 
in detail e.g. in the renormalization scale in a dimensional scheme 



(e.g. MS) serves to separate short- distance from long-distance physics. The 



Appelquist-Carazzone decoupling theorem can be implemented prop- 
erly in the (mass-independent) MS scheme by hand by matching the effec- 
tive theories below and above thresholds. The advantage of having a mass- 
independent scheme is that the renormalization group /3-functions do not 
explicitly depend on the scale /x, while the validity of the decoupling the- 
orem guarantees that all intuitive reasoning based on a so-called physical 
renormalization scheme remains still true. 

When the effective theory contains two heavy mass scales mi, 1712 of compa- 
rable magnitude, it is usually a good approximation to integrate out both 
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particles at a common scale. On the other side, if the ratio x = (mi/m2)^ is 
very small (i.e. x <^ 1), even if the coupling constant is small, the product 
a In X may become of order unity, and one is then forced to sum all powers of 
this product, while corrections to the sum are suppressed by powers of a or 
X. Sometimes the situation is less favorable and lies somewhere in between, 
as is the case e.g. for the SM with a heavy top quark |jl3[ of, say, 175 GeV. 



For the process b —> 57, the most important correction is the QCD running 
between the W and the b scale, whose size is (parametrically) given by 

asimw) ■ ln{m^^/mbY ^ 0.7, 

while 

a5'(mw) ■ ln(mt/mw)^ — {m-^/mtY ~ 0.2 

indicates that one might miss numerically important pieces if either of the lat- 
ter would be neglected; compared to next-to-leading order corrections which 
are of order as — 0.1. i What one can achieve with reasonable effort is to take 
into account the resummation of the leading terms in the limit of a heavy top 
quark, and then simply adding in the nonleading terms, thus neglecting terms 
which are (up to logarithms) 0{a{mY^/mtY). The choice of scale where to 
add these nonleading terms is at this stage completely arbitrary and can 
only be answered by a calculation of the power corrections. The remaining 
uncertainty is, however, less important than neglected next-to-leading order 
corrections. 

Let us now turn to the application to the 6 — > 57 transition. The effective 
Hamiltonian of interest may be written as a sum of Ai? = 1, AS = 1 
operators: 

A suitable operator basis {Oi} will be given below. 

In general the definition of the operators in (|l|) will require the specification 
of a renormalization scheme. From the fact that the effective Hamiltonian is 
independent of the renormalization scale, one derives renormalization group 



•^Of course a full calculation of the next-to-leading order corrections is necessary to 
resolve the well-known ambiguity in the choice of scales in leading-order calculations. 
This would require however the computation of three-loop anomalous dimensions for the 
process under consideration. 



4 



equations for the composite operators Oi{^) and the coefficient functions 
Ci{ii). The renormahzation of a composite operator is formally defined in 
terms of the divergent renormahzation constants Zij which relate renormal- 
ized and bare operators: 



(2) 



Since the bare operators are /i-independent, the renormalized operators de- 
pend on the subtraction scale via the fi dependence of the Zij: 



djjL 



O 



bare 



-likOk, 



d 



where 

lik 

is the so called anomalous dimension matrix. 



(3) 



(4) 



From the scale independence of the effective Hamiltonian (|l|) one derives the 
renormahzation group equations for the Wilson coefficients C,: 



E(7n,^.(.) 



(5) 



If QCD corrections are neglected, the solution to this differential equation 
is straightforward. When QCD corrections are included, it turns out to be 
favorable to eliminate the derivative with respect to the renormahzation scale 
in favor of a derivative with respect to the coupling constant: 



(6) 



Here (and in the following) denotes the QCD coupling constant, and (3 
fi^dg^/dfi) is the QCD beta function. 

The solution to the differential equation (P) is then given by 



93 (m) 



exp 



dg 



93 (W)) 



Pig) 



(7) 
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where Tg means an ordering in the couphng such that g increases from right 
to left (for /i < Ho). Since our anomalous dimension matrices will be times 
a purely numerical matrix, 

7 = ^^7 + 0(^3^) 
the (yf-ordering is superfluous, and the (yf-integration is trivial: 



C(/i) = exp 



2b asifio) J 



Cipio) (8) 



The most convenient way to calculate the anomalous dimension matrix 7 is 
to consider Green functions with insertions of composite operators. Denote 
by r^,^ a renormalized n-point IPI Green function with one insertion of the 
operator Oi. The anomalous dimension jij that determines the mixing of Oi 
into Oj may then be simply read off from the renormalization group equation 

for rg, 

7.rg = - (/^l; + + Imm^ - n7e.,) T^g (9) 

Here 7^ = {^/m){dm/ d^), and n7ext accounts for the wave-function anoma- 
lous dimensions arising from radiative corrections to the external lines of the 
Green functions. 

We shall use dimensional regularization with minimal subtraction (MS), d = 
4 — 2e. The SU(3)c x U(l)em covariant derivative then reads: 

D^ = d^- ^/i^^^sG^X" - tfi'eA^Q (10) 



We will use the background field gauge throughout this work. The 



anomalous dimensions of the fields are in this case given by: 

-I A -AA -I 

7quark „ o 9 ' 7squark o o 9 ' 7gluon ; 

3 Svr^ 3 8tt^ gs 

7m = -4f^, 7^ = -6^, /3 = 6f^. (11) 
St:^ Stt^ Stt"^ 

The coefficient appearing in the /3-function has the value 

11 1 1 
^= -y + 3^/ + Y^^-g + ^-s, (12) 
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where nf,nq and Ug are the number of active quark flavors, squarks and 
gluinos, respectively. 

In all cases considered below, the operator basis of choice will contain the 
following set of operators involving only light degrees of freedom (i.e. photons, 
gluons, and "light" quarks with masses below mw):^ 



Dimension d + 1: 



Dimension d + 2: 



0} 



LR 



>1,A 



Pi 



167r2 

^3 



i 



167r2 

eQh 



IGtt' 



IGvr 



The tensors appearing in Pl'^, A = 1, ... 4, are defined by: 



(13) 



rri3 

jxua 



T. 



fiucr 



T. 



/ii/aT 



7 75 



(14) 



In order to apply the procedure outlined above to the MSSM case, we will 
consider in a first step the extension of the calculation by Cho and Grinstein 



r^l to the case of a type-II two-Higgs doublet model. There are already two 



cases to consider, namely that the charged Higgs can be either much lighter 
or much heavier than the t quark. We shall then explain in detail how one 
adds to this picture the contributions induced by the chargino loops. 



^Note that our normalization differs from ref. We have omitted their operator P^, 
since the corresponding Wilson coefficients will always be zero, and none of the operators 
under consideration will mix back into it. 
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3 Two-Higgs doublet model 



3.1 rrit > mH± 

If the top quark is heavier than the W and the charged Higgs, then the first 
step is to integrate out the top quark at the scale ^ = m^. This leads to 
an effective field theory for /i < mt without the t, but with new vertices of 
dimension larger than four that contain the virtual t effects. For the process 
under consideration, we need, in addition to the operator basis (p!3D, further 
operators. 

In general, in the range mt > /i > mw, one has to consider higher dimen- 
sional operators which contain the W's, the would-be Goldstone bosons (j)± 
and the charged Higgs field. By naive dimensional analysis, we expect that 
higher dimensional operators are suppressed by inverse powers of the ratio 
Xtw = {'fntlrriY^Y . Since this ratio is not very large for phenomenologically 
acceptable top quark masses, the effects of the higher dimensional operators 
are not necessarily small, compared to the leading dimension 5 and dimen- 
sion 6 operators. Also, the matching conditions at threshold in general are 
combinations of rational functions and polynomials in Xtw- 

Nevertheless, we shall take the approach motivated in the previous chapter 
and keep only the leading operators and the leading terms in the matching 
contributions. Although we are unable to calculate the power corrections, 
we shall later add the sub leading terms in 1/xtW; so that we get the same 
result in the limit of neglecting strong corrections for /i > m-w, as when all 
heavy particles are integrated out simultaneously at the W scale. 

The relevant part of the interaction Lagrangian in the charged current sector 
reads 



where U = {u, c, t) and D = {d, s, h) represent up-type and down-type quarks, 
respectively, Mu = diag(mu, mc, mt), Md = diag(md, mg, mb) are the quark 
mass matrices, g2 = e/sin6'w is the gauge coupling of SU(2)w, -Pl,r are 
projectors on the left- and right-handed components of the fermions, and K 




+ 



(15) 
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is the Kobayashi-Maskawa matrix. In the present work, we shall neglect the 
masses of the quarks of the first two generations whenever appropriate. 

From these expression one can see that the leading terms for Xtw ^ 1 come 
from vertices which involve the charged would-be Goldstone bosons 0± and 
the top quark, since they are proportional to the top quark mass. For this 
reason, in the range rrit > fi > itlt^, we shall need, analogous to the findings 
in 1 13], the following operators with external would-be Goldstone bosons, in 



addition to the operator basis ( [TB| ) 



QhR 


= u — 


Rl 


= ^^'^ 

^ 167r2 


Rl 


■ 2e 93 
= Ifl ^ 


Rl 





Sl7m^R (16) 



The inclusion of explicit factors into these operators is motivated by the 
Oilman- Wise trick [^, which allows us to have all one- loop contributions to 



the anomalous dimension matrices to be of 0{gf), so that the diagonalization 
of these matrices is scale independent. We will freely use this trick later on. 

The interaction Lagrangian for the charged Higgs with the quarks reads: 

Lh±./, = (h+U [cot (3MuKP^ + tan PKMdPr] D + h.c.) , (17) 

V2mw 

with tan/3 = Vi/v2 being the ratio of the vacuum expectation values of the 
Higgs fields which give rise to the masses of up- and down-type quarks, 
respectively. 

The interaction ( [T7| ) has the same structure and quark mass dependence of 
the couplings as the interaction of the would-be Goldstone bosons 0±, see 
(|T5|). In the limit XtH = ('^t/'^H±)^ ^ 1, keeping only the leading terms in 
1/xtK, we are lead to the following operators with charged Higgs bosons we 
have to add to our operator basis in the range nit > fi > "^h* • 

iDvr^ 
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= Z/i^^ 


9l 


Rl' 


= Z/i^^ 


9l 


Rl' 




9l 
167r2 



{D^H^)H- (18) 
H+{D^H-)s^-f,bn 



3.1.1 Matching at fi = rrit 



For fi > rrit, our effective theory is a fully renormalizable theory, which still 
contains all particles and interactions, so in this case all coefficients of our 
effective Hamiltonian are zero: 

Ci(/i = m+) = for all i (19) 

When we cross the t threshold from above, i.e. when we integrate out the 
top quark at fi = rrit, we obtain the following changes to the coefficients 
of the effective Hamiltonian, due the interactions from the would-be Gold- 
stone bosons from matching the three-point functions F''^'^ and F^'^^ [|l3l (here 

1 



2 



11 

18 



Q 



4gb 

AC- 
^{0,iy)<--QLR - — -0- 

93 

^{4>,w)Cr3 = (20) 
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Similarly, there are contributions from the interactions with the charged 
Higgs bosons:! 



1 

2 



Qt 



1 

2 



11 



A 



(H)Cp2 



A(H)Cpi' 



A(^^)C^3' 



— cot^/J 
9 ^ 

-QbA(/f)Cp4 

3 



cot^ /? 



4Qb 
167r2 



cot^/3 



A 




:5— cot p 

9i 



(21) 



At this point it is worthwhile to note that, had we not matched at the scale 
= rrit but at a different scale (or used a different subtraction scheme), we 
would have found logarithmic contributions in the matching corrections to 
the coefficient of P^- 



A(<^,Ty)Cp2 = — 



4 6 



m; 



A 



3^1, /i 
4 + 6^"^ 



cot^/3 



(22) 

These logarithms which vanish for /i = mt are regenerated at lower scales by 
the renormalization group for the effective theory below mt. It is therefore 
not surprising that they are present in the full expressions for this coefficient 
given in the appendix, when both particles in a loop are integrated out at 
the same scale; there it appears as an unsuppressed logarithm of the mass 
ratio of the particles in the loop. 



^We prefer to keep the contributions from each interaction separately, for there are 
different cases below. 
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3.1.2 Running below rrit 



The anomalous dimension matrices for the mixing of the operators Oi and 
Pi has already been given in ref. [|13i- For completeness, we quote the result 



obtained in this work. 

First, there is a mixing of the operators Q, R with would-be Goldstone boson 
fields into the operators without (O, P): 



(23) 



Note that this mixing back is of order 0{g^) due to our choice of the co- 
efficients in front of the operators Q, R, and not due to "proper" QCD 
corrections. 

For the QCD-induced entries in the anomalous dimension matrix, one has 





Olr 




Pi 


Pi 


Qlr 


( 








\ 
















Rl 








l/6Qb 





Rl 


I 





-i/egb 


/ 



7 



01k 

01k 

>1,1 

L 

)1,2 
L 

,1,3 
L 

Pi 
Pt 



PI 
Pi 
P' 



OIk 

I 20 

3 





6 
4 
2 



V 



•-^LR 

1 

2 
3 


2 

3 
2 

1 

1 

2 






3 
16 





1 

2 


4 
3 



Pi:' 


pl,2 
^L 


pl,3 
^L 


pl,4 
^L 






































2 

3 


2 


-2 


-2 


113 


137 


113 


4 


36 


18 


36 


3 


-2 


2 


2 
3 


-2 


113 


89 


113 


4 


36 


18 


36 


3 



























PI 






9 
4 



9 
4 






PI 

\ 








} 



(24) 



Similarly, the mixing among the operators with would-be Goldstone boson 
fields is known to be: 



7 



Rl 
Rl 
Rl 



QhR 


Rl 


Rl 


Rl 


-2b 













-2b 














-2b 














-2b) 



(25) 
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Obviously, the same mixing matrices are found when one considers the mixing 
of the operators with charged Higgs fields, i.e. when one replaces Qlr Q'lr, 
Ri Rl' in eqs. (|2|) and (^. 



3.1.3 Matching at /i = mH± and n = mw 



In the process of scaling down, when we encounter the charged Higgs or 
W threshold, we have to integrate out the or W and would-be Gold- 
stone bosons, respectively. Due to decoupling that has to take place below 
threshold, we shall remove the operators Q', R' from our operator basis for 
/i < and Q, R for /i < mw Again we obtain the finite changes of the 
coefficients of the operators O and P by matching Green functions calculated 
in the theories above and below threshold. 

Since we neglect small terms proportional to ttt-u or ttt-c, we find no nonva- 
nishing contribution from the matching of the effective theories above and 
below n = mH±, i.e. our Wilson coefficients are continuous: 

=C,(m-±) (26) 



Matching the effective theories above and below /i = m-yv, we find the fol- 
lowing changes in the coefficients of the effective Hamiltonian (here AC = 



((l>,w)Co'f 













= A 


^{<l>,w)Cpi,i 




2 

~ 9 


^{<I,,W)C pl.2 


7 

~9 




^{<l,,W)CpiA 


= 1 






1 






3 





(27) 

Again, had we matched at a different scale /i ^ Tn-^, we would have found 
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different matching contributions for some of the coefficients: 





2 , /x^ 



9 3 



A(0,vy)Cpi,2 



1 A , /x^ 



) 



(28) 



But the coefficients of In/x^ are just the coefficients of those logarithms in 
eqs. ( [77| ) which give the leading (divergent) contribution to the Ci in the 
limit of small quark masses. These logarithms are regenerated by the renor- 
malization group running in the low energy effective theory valid at scales 
yU < mw and therefore need not be discussed here any further. 

We shall to now use our freedom to add subleading terms in 1/xtW; l/a^tH 
to the coefficients Ci. In order to see how this is accomplished, let us for 
the moment neglect the proper QCD corrections, so we have to consider only 
the entries in the anomalous dimension matrix given in (^B]). Solving the 
renormalization group equations (|^), we find that only one coefficient runs 
below mt. 



where the first term in parentheses is due to the mixing of into Pl, and 
the second due to . We see that the renormalization group reproduces 
the logarithmic terms already discussed in eq. (p2|) , which would have been 
there, had we done the matching at a different scale. 

The subleading contributions are found by taking the standard one-loop re- 
sult from integrating out both particles in the loop at the same scale, see 
appendix, and subtracting the leading contributions that we have found from 
the matching contributions (pO| , PD and the running ( ^91) without QCD. We 
will always refer to this procedure for obtaining the subleading terms in the 
rest of the present work. 

Let us for the moment assume that > mw When we integrate out 

the charged Higgs at /i = mH± , we obtain the subleading contributions from 




(29) 
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(gT]®|Il|): 



= (2F2(xtH) + + 2F4(xtH)) - ^) cot2/3 

/ 2t S \ 



V 3 ^ ^ tn; my 4V my; g 

A'(H)CpM = -QbA'(^)Cp4 = (xtHF4(a;tH)-^)cot2/3 (30) 



A;„A, = i-(..„(lf,(..„) + ^.,(..„))-?-_J^)cot=;3 

The functions Fj(x) are given in appendix A. One may easily verify that 
the terms on the r.h.s. are of order O^l/xtn), so they are truly subleading. 
Especially there is no (leading) logarithmic dependence of the matching con- 
tributions to Cp2 on the mass ratio Xtu, since all such dependencies must 
come from the renormalization group. 

Of course there is an ambiguity in the choice of scale where to add the 
subleading contributions. This ambiguity can only resolved by computing 
the power corrections, which fortunately differ from our naively adding the 
subleading terms only by a next-to-leading contribution. We shall define our 
procedure by assuming that taking the scale equal to the mass of the lightest 
particle in the loop is a suitable choice. 

After scaling down from /i = mH± and adding in the leading matching con- 
tributions at /i = mw, we will also consider the subleading contributions. In 
analogy to the previous case we find: 



2 
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= ^^^^(2F2(xtw) + i^3(a:tw) + 2F4(xtw))-^ 
^U,w)Cpi,^ = ^^^^^^(F2(xtw)-F3(xtw)-2F4(xtw)) + ^ 



In Xtw 



^[ct>,w)^Pl = ^ ((^tw + 2) QF3(xtw) + i^4(a;tw)) _ 

^y)Cp4 = 2xtH-^3(a;tH) - 5xtH-^4(a;tH)^ (31) 



3.1.4 Reduction by equations of motion 

In order to be able to use the results from previous calculations for the 
running between the W and the b scale, we have to match our operator basis 
to the operator basis employed there. To this end, we use the equations of 
motions, as in ref. |TB[. For the effective Hamiltonian just below the W scale, 
one then finds: 



eff 



AC 

-^K,^K E Qim^mm^) (32) 



EOM 4G'f j-^ 



V2 



Cnl + CnlA Cd1.2 H tCtjIA ] Ot R 

-Col + C*o3 — -Cpl.l — -Cpl.2 + -Cpl,4 — -Cp4 
2 ^LR ^LR 2 L 4 L 4 L 4 L 



H — (four-fermion operators) 

1677^ 

Since we are only interested in the leading contributions from the QCD cor- 
rections due to a large mass splitting, we may drop the contributions to the 
four-fermion operators in (|32D since these are suppressed by a factor gl/lQ-n'^ 
and therefore truly nonleading. 

To this expression we have to add of course the standard four-fermion op- 
erators (&L7^SL)(gL7^gL), q = u,c (with the appropriate CKM mixing coef- 
ficients) from integrating out the W. The Wilson coefficients obtained this 
way at the W scale may then be used as input for the renormalization group 
running R P, IT2[ down to the b scale. 
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3.2 mH± > TTit 



If the charged Higgs is heavier than the top quark, the picture becomes a 
httle more involved. As we run down from large scales, we first encounter the 
threshold of the charged Higgs. Therefore, as a first step, we integrate out 
the charged Higgs. In the same spirit as in the previous case, we shall now 
be mainly concerned with the leading contributions in the limit Xtn <^ 1- 

In the range mH± > fj, > rrit, after integrating out the charged Higgs, we 
have to deal with four-fermion operators of dimension 6 that involve a b, an 



s, plus a quark-ant i- quark pair. Besides the operators (|l^), our operator 
basis contains: 



^1 = 

^2 = 
^3 = 


q 


^4 = 




^5 = 




^6 = 





loTT^ rrit 

Sio = i/^'^lil^(^"^.^^R)(^t^'^^&^) (33) 

Here a,f3 are color indices of the quarks, and the sums run over all active 
flavors. Again, the inclusion of the additional factors is motivated by the 
Gilman-Wise trick [^, as are the factors mh/mt to keep the anomalous di- 
mension matrices mass independent. The different normalization of 5*1 ... 5*6 
and 5*7 .. . 5*10 will be explained below. 
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Integrating out the charged Higgs at /i = mH±, we find to leading order in 
Xtn = (mt/mH±)^: 



Cq^ 



Cpl,4 



C p2 



1 



Cpi,3 = XtH ■ cot^ /3 
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9 

XtH ■ cot^ /3 



-XtH ■ cot^ /3 



i-(-ix,H)cot^/3 

--XtH cot^/3 

167r2 
— o~3;tH 

9i 



0, 



1,3... 7, 9, 10 



(34) 



Let us start again with the mixing back of the operators S into the operators 
O and P. Because of the chirahty structure of the operators, we find two 
different situations at one loop. The operators 6*1, ... appear to have a 
zeroth order mixing ((73) at one loop into the operators P. 
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(35) 



However, by inspection of the equations of motion (|32D one sees that the back 
mixing vanishes to this order; therefore we may simply drop this contribution. 



18 



As is well known, one has to consider this mixing at two-loop order. The 
anomalous dimension matrix can be derived from eq. (25) of ref. Q, and 
reads in our normalization 



r)2 



7 



Si 
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3 
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S2 
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54 


S3 




i + |n, 


S4 






s. 




-1 - ¥f 


Sq 
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54 'v 





224 

27 
232 

27 

27^/ + 4% 
16 



Here n 



f 



27^/ - J 

+ n-d is the number of active flavors, and n f 



(36) 



nd + {Qu/Qd)nu. 



On the other hand, the mixing of 5*7 
vanish at one loop: 
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. . S'lo into the operators O does not 
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IQt, 
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(37) 



Again one may verify that these entries of the ADM are consistent with the 



In /i-dependence of the matching contributions (@). 

Let us now turn to the mixing among the four-fermion operators. Since 
the considered operators are all of dimension d + 2, and because the QCD 
interactions preserve chirality, the operators S'i^...^6 and the operators S'7^...^io 
will mix only among themselves, respectively. 

The one- loop mixing among the S*! ... 5*6 is well known ||3^ : 
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(38) 
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For the mixing of 6*7 . 



S'lo "we find: 
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(39) 



Since the operators Olr are dimension c? + 1, there is no mixing back into 
'S'i,...,io- 



Note that with our chosen normahzation of the operators as given in ( p3|) all 
relevant mixing occurs at order g"^, and all entries in the anomalous dimension 
matrix are dimensionless. 

After running down to /x = mt, we integrate out the t quark. As far as the 
operators Si, S2, 5'7,...,io are concerned, they are just removed, since they give 
no contribution to the matching; for the operators S'3,...,6 the t quark has to 
be excluded from the sum, because it is "inactive" for n < mt. Again we 
will take into account the sub leading terms in x^h according to the general 
prescription given in section 3.1.3| . Then we will continue as in the case for 
the Standard Model with a heavy top, except that the coefficients Ci{mt) 
are now nonvanishing. 



4 Super symmetric contributions 

4.1 Flavor changing chargino interactions 



Let us in analogy to JI^ denote by W^, iff and H2 the superpartners of 
the W and the charged components of the Higgs fields, respectively. Define 
the two component spinors ^/^^ by 



1,2 



(40) 



The mass term for the W-inos and higgsinos then takes the following form: 

Lm = -^~Mtl)^ + h.c. (41) 
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where the mass matrix is given by 

j^^f v^mwsin/3y 
Vv2mwcos/5 fJ-h J 

where M2 is the soft SUSY breaking mass term for the W-inos at the weak 
scale, and /ih is the renormahzed Higgs mixing parameter. 

This mass matrix may be diagonahzed with the help of two unitary matrices 
U, V such that 

U*MV^ = M^ = diag(mi, ma) (43) 

is a diagonal matrix with nonnegative entries. The corresponding charged 
mass-eigenstate 4-spinors are the charginos 

We shall find it however more convenient to rewrite the interactions of the 
charginos by their charge conjugates 

x;^{xty=c(xtr={^^), (45) 

so whenever we refer to charginos below, we mean the x7 given in (|45|). 

Let us apply these definitions to the interactions of the charged gauginos and 
higgsinos and convert to 4-spinor notation. Neglecting for the moment the 
mixing of quarks and of squarks and concentrating on the terms involving b 
quarks, the relevant Lagrangian for char gino- quark- squark interactions reads: 

hbt = -92V:JUx^PLb) + 92XtV^il{x^PLb) + g2KUdi{x^PKh) + h.c. , (46) 
and the couplings Ag are proportional to the Yukawa couplings: 

= ^ ~^ ' ^''^^ R- 

V 2mw sm p v 2mw cos p 

A similar expression is found for the interactions of the charginos with the 
quarks and squarks of the second family. In this case one can however ne- 
glect the terms proportional to Ac, As, which originate in the coupling of the 
higgsino components of the charginos to the quark and squark fields. 
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Since the Yukawa couplings of the matter fields to the Higgs fields are not 
fiavor diagonal in a weak interaction basis, we have to take into account the 
mixing among quarks and among squarks. Let us denote as in ref. |]T^ by 
qiL,R the squark current eigenstates (where q = u,d, and / = 1,2,3 is the 
generation label), (ja {a = 1, . . . , 6) the corresponding mass eigenstates with 
masses rha- We define the 6x3 squark mixing matrices Fql r by 

Ql,r = T^ql,kQ ■ (48) 

The relevant chargino interactions involving down-type quarks may then be 
written as 

L^du = -92 E K X, {G^'^'Pl - H^^'Pr) di] + h.c. , (49) 

where 

= v;,T^^^-V*,iTuKAuKf 
W^' = U.^iVui^kDY' (50) 

Here Kjj = Mf//(^/2mwsin/?), = Mr,/ {\/2my^ cos (3), are proportional 
to the Yukawa coupling matrices for up- and down-type quarks, respectively. 
Note that we neglect the masses of the light quarks, and therefore we set the 
Yukawa couplings of the light quarks to zero. 

Since we are interested in the b ^ s'j transition, we find it convenient to 
define 

g^-^ = -— for I = b,s- W-^ = -— (51) 
J^ti J^th 

Unitarity of the mixing ( ^5]) implies that 

6 6 

E-pai p*afc r \ ^ -pai -p*ak pi fKO\ 

^ QL,R^ QL,R — (Jik , 2^ J- QL,R^ QR,L — ^ ■ \^^) 
a=l a=l 

and therefore 

6 6 
a=l a=l 



22 



After having described our conventions, let us now turn to the evaluation of 
the QCD correction. As the squarks and the charginos can have large mass 
splittings, the procedure of matching and running will become more involved 
but still remains straightforward. We will give all ingredients, but the precise 
procedure will depend on the details of the spectrum. 



4.2 Effective operators from heavy squarks 

In the process of evolving down, if we encounter the threshold of an up-type 
squark Ua, we will integrate it out. This generates effective four-fermion 
operators made out of the quarks b, s, and the active charginos Xj- We 
extend our operator basis by the following operators (no sum over j): 



Wl = /x^^^(sl7^&l)(x^r7.x{.) (54) 
The matching contributions at /i — rha for rhj ^ rha are: 

A(.„)Co3 = g*.^n^^^.^f-^^V (-1) 
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= ^--^-'^ • (^) '•(-!) (55) 
The mixing back of these operators into the O's and P's is found to be: 

Oln ^l'^ pi Pi 

\ 

-l/Qb (56) 

-2/3Qb ; 

Since the charginos carry no color charge, the renormahzation of these oper- 
ators is particularly simple, 









wl 




( -2b 
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7 = Wl^ 





f -26 





Wi 


I 





-2b J 



and there is no mixing of the O and P operators back into these. 

If we cross the threshold of chargino Xj at = rhj, the operators will just 
be removed; they do not give any matching contribution to leading order. 



7 



wli 

Wl 










4.3 Effective operators from heavy charginos 

Let us now consider the case that we encounter the threshold of chargino Xj 
at yU = rhj. If there are still active up- type squarks Ua, we have to extend 
our operator basis by the 2-quark-2-squark operators (no sum over a): 

Dimension d+ 1: 

QIr = m^uMsZb'^ 

Dimension d+ 2: 
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^l'" = mftif Tr 

Rt^ = i{D^Uaru^,sl-f,hl (58) 

For a running over each active up- type squarks and using Zja = {rhj/rha)'^, 
we find the following leading matching contributions (i.e. for large Zja) at 
H = rhf 



L 



rhj J V 9 / 



'mw\^ / 3\ 1 



2 



g*jas^jah . ^ / . /gN 

mb V "^j / 



^x,)% = A(^.)C^n,. = 71 = 3,4,5,6 (59) 
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A straightforward calculation gives for the back-mixing at one-loop order 
(but order in our chosen normalization) 
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(60) 



Qlr 
Qlr 

Rl" 

Rf 
Rl" 



Again one sees that that, similarly to the case of the four-quark operators, 
the mixing into the magnetic moment operators vanishes after applying the 
equations of motion. Therefore we have to consider this mixing at two-loop 
order. 

The actual two-loop calculation of the mixing of the two- quark-two- squark 
operators ( |58D into the magnetic moment operators is performed analogously 
to the corresponding calculation with insertions of four-quark operators (see 
e.g. 0). In figure Q we show the relevant diagrams and one-loop counterterms 
contributing to the mixing of the two-quark-two-squark operators into the 
operator Olr. As we prefer to work off-shell, we have to consider only 1- 
PI diagrams. The main advantage is a simplification of the extraction of 
the divergent parts of interest by focussing on the coefficients of the tensor 
structures that are defined by our basis (0). 

Similarly to the corresponding calculations with insertions of four-fermion 
operators, using the equations of motion (^) greatly reduces the compu- 
tational effort. Figure ^ shows typical diagrams which do not contribute 
because their sum can be shown to be proportional to (7^?^ — 5'/^^), and 
therefore need not be calculated. 

In the case of the mixing into Olr, due to the non-Abelian interactions of 
the gluons, we have to consider the additional diagrams and counterterms 
shown in figure ^. 



26 




Figure 1: Feynman diagrams contributing to the mixing of the 
two- quark-two- squark operators (|58| ) into the operator Of^^. A full square 
denotes the insertion of a two-quark-two-squark operator, while an open 
square denotes a one-loop counterterm. Diagrams which are related to the 
ones above by reflection are not shown. 



We obtained the following mixing coefficients (A^ = 3): 
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In addition we need the mixing among the two-quark-two-squark operators, 
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Figure 2: Feynman diagrams whose contribution to the mixing vanishes after 
apphcation of the equations of motion. 



where the squarks are of the same kind, 
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(62) 
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Figure 3: Additional Feynman diagrams which contribute to the mixing of 
the two-quark-two-squark operators into the operator Olr. 



and for different types of squarks {a ^ b): 

\] (63) 

36 / 

Note that there is a mixing of some of the two-quark-two-squark operators 
into four-fermion operators: 

(64) 

If there are squarks hghter than the top quark, we also have to take into 
account the mixing of the four-fermion operators into the operators R: 
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In principle there is also a QCD-induced mixing into operators with two 
quarks and two down-type squarks, which we also would have to include if 
we were considering the contributions induced by gluinos and neutralinos. 
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In most scenarios, the mass splitting of down-type squarks is much smaller 
than for up-type squarks. For the supersymmetric contributions to be nu- 
merically relevant the lightest squark (which is usually the lightest stop) 
must be significantly lighter than the other squarks. As has been argued in 
the introduction, contributions from these operators are strongly suppressed, 
and also inclusion of these operators into the mixing would lead to only a mi- 
nor effect compared to other neglected corrections. Furthermore, all Wilson 



coefficients that contribute to mixing via (|63| , |64D are proportional to cot^ f3 
and therefore suppressed in the large-tan /? limit. 

Again, if we cross the threshold fi = rfia of squark Uai the matching con- 
tribution vanishes to leading order, and the operators O*^, are simply 
removed. We will also add in the corresponding subleading contributions 
each time a pair (a, j) of squarks and charginos has been integrated out, i.e. 
aX, II = min(ma,mj). 



5 Results and Discussions 

As the full anomalous dimension matrix is quite large and changes its struc- 
ture every time we cross a threshold, it would be a big effort to diagonalize it 
in every step. It is much simpler to directly evaluate the solution of the 
RGE numerically. Before we proceed, let us comment on some simplifications 
that result from the use of the equations of motion, since we are eventually 
only interested in the coefficient of the magnetic moment operators at the h 
scale. 

First we note that the operators Qlr and R]^, and in the case of mH± < 
rrit their "primed" counterparts, which appear in intermediate stages of the 
calculations, turn out to be superfluous as they do not give any contribution 
in the process of matching, nor do they mix into any other operator. Second, 
although the coefficient of does get matching contributions, and many 
operators mix into it, it can be ignored, since it vanishes after applying the 
equations of motion. Third, the operators i?L and i?L only into Pl, 
which vanishes by equations of motion, and may therefore be omitted from 
the beginning. Extending this reasoning to , R'l' , W^lr? and Ftf^"' 
shows that they may also be disregarded. 
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Next, one may convince oneself that the apparent zeroth-order mixing of 
some operators (see eqs. (|35D , (|60D ) vanishes after apphcation of the equations 
of motion, so all mixing occurs at order (^fl/Svr^), as promised. 

Let us first rediscuss the effect of the QCD corrections to the Standard Model 
contribution. For the contribution from the W-t-loop there is a QCD en- 
hancement of the coefficients C?J1 (mw) and (mw) of the order of 10- 
18% and 15-22% for rrit = 130 . . . 250 GeV, respectively, which after scaling 
down to fi = rrih and including the contribution from the four-fermion op- 
erators leads to an additional enhancement of the decay rate within the SM 
of the order of 12-23% ||13[, compared to the case when both t and W are 
integrated out at /x = mw- This large correction, which seems to compare 
quite well with the naive estimate given in the introduction, is a confirmation 
that a full next-to- leading order calculation is quite important. 

The magnitude of this effect may be understood by solving the renormaliza- 
tion group equation for the leading terms. After application of the equations 
of motion, their contribution turns out to be quite simple: 



^oi(-w).,, = f:^#^V'''' + ^ + (subleading) 



SM 24 \a3(mw) / 3 



Co^ (^w) = ^ ( ^^T^ 1 - ^ + (subleading) (66) 



LR 



14/23 



SM 3 \ a3(mw) / 3 



The reason that (^) leads to positive corrections is essentially due to the 
fact that the effective matching contributions at /i = mt (pO[) and at fi = mw 
( ^TD have opposite sign (which is a remnant of the GIM mechanism), and 
therefore lead to coefficients of opposite sign but comparable magnitude of 
the first two terms on the right-hand sides of (|66D . It has long been known 



3^ that the QCD corrections tend to soften the GIM-cancellations between 
different up-type quarks if they are nearly degenerate; but even for a heavy 
top quark (i.e. mt ^ m-w) there remains a finite enhancement, as can be 
explicitly seen from these expressions. Note that (|66D gives only the leading 
terms; the subleading terms being suppressed by only a factor of (mw/m-t)^. 

Next let us turn to the contribution from the loop with a charged Higgs. 
For this case we have solved the renormalization group equation numerically. 
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Figure 4: Correction ( |67D to the coefficients of ^^^(mb) from a loop with 
t quark and charged Higgs with leading QCD corrections from large mass 
splitting to the case when both particles are integrated out at /i = mw The 
mass of the t quark is assumed to be 175 GeV. The dotted, long-dashed, 
dashed and sohd line correspond to tan/3 = 1.5, 2, 3, and 10, respectively. 

using as input parameters: 

mb = 4.5 GeV, mt = 175 GeV, 
mw = 80.22 GeV, a3(mz) = 0.123 

The resulting correction 

5h± 



^eflf, naive / \ 



- 1 



(67) 
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to the "naive" result, obtained by integrating out t and H='= simultaneously 
at the W scale, is shown in figure H for mw < fnii± < 750 GeV and tan /3 = 
1.5, 2, 3, and 10. At sufficiently large tan/? (i.e. tan/? > 3), the correction 
turns out to be essentially independent of tan j3. This is quite understandable 
since the tan /^-dependent pieces are actually proportional to cot^ p. 

For a light charged Higgs, i.e. mH± < rrit, there appears to be a further 
reduction of this contribution compared to the naive result. Indeed, in the 
limit of large tan /?, and assuming there is no light squark or gluino with mass 
below mt, one finds the following simple analytical result for the charged 
Higgs contribution, valid for rri]^ < fi < mt. 



contrary, the leading coefficients get suppressed as they are run down from 
the mt, compared to the subleading terms which (according to our discussion 
in section |3.1.3|) get only suppressed by the evolution from /i = mH± down 
to fi = mb. For sufficiently small mH±, the additional QCD corrections are 
then essentially due to the running from /i = mt to /i = mH± . Note that our 
"corrections" are counted relative to the case when both particles in the loop 
are integrated out at the common scale fi = m-w, which is obtained from 
( p5D by substituting a;3(mt) — > as^niy^). Thus, for mH± < mt, integrating 
out t and at the t scale appears to give a more accurate result than at 
fi = mH± or fi = mw- 

On the other hand, for mH± > mt we found only a minor suppression of a 
few percent, which is essentially the result of a partial cancellation of the en- 
hancement coming from the scaling between mH± and mt (due to one negative 
eigenvalue of the submatrix (|39| ) for the mixing of four-fermion operators), 
and of a reduction from the scaling between mt to mw Unfortunately, we 
were unable to obtain a simple analytical solution for this case. 

In the case of the chargino contribution, things are more complicated, since 
one has to consider in general the dependence of the amplitude as a function 
of several parameters, namely the mass spectrum and the mixing angles 
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for the charginos and the up-type squarks. However, it turns out that the 
essential features may already be studied for the case of sufficiently large 

m. 1231, m 



In 



tan/3, which is in the center of recent interest 
this case, the parameter Ab (^Tj) may become of the same order of magnitude 
as the parameter At- Assuming furthermore that the mixing in the squark 
sector is essentially the same as in the quark sector, which is quite natural 
in supergravity models where the soft SUSY-breaking is characterized by a 
common scalar mass at some unification scale, the quantities Q and Ti, as 
defined in (^l]), are then necessarily of the same order of magnitude, the 
terms proportional to the ratio rhj/nih will dominate the amplitude, and the 
corrections become tan /3- independent. 

In this particular limit, one can find an analytical result for the leading 
terms. For the case of the chargino being much lighter than the squark, the 
coefficients read: 
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while for the other case of a squark much lighter than a chargino, we get: 
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(70) 



where now 



m 



w 



2b yasirhj) 



At first sight the terms proportional to l/a^ in (|69| ) might be embarrassing, 
but a closer look shows that their difference is (to leading order) just some 
number times Inlfi/rha) and therefore finite in the limit — > 0. 
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Unfortunately, the interpretation of these expressions is aggravated in both 
hmiting cases since the number of free parameters in the general model is 
quite large, and due to eqs. ( p^D one has a supersymmetric version of the GIM 
mechanism, which leads to a partial cancellation of the leading terms under 
consideration. This renders it difficult to estimate the actual corrections due 
to the mass splitting between charginos and squarks by using (|K 



or (|7g). 

Some features of these expressions may still be studied under the following 
assumptions: i) the squarks of the first two generations are degenerate with 
mass ifiui h) the mixing in the squark sector is the same as in the quark 
sector (i.e. the gluino- quark- squark couplings are fiavor- diagonal even in the 
mass eigenstate basis), and iii) the mass matrix for the stop is given (in the 
(tL,iR) basis) by 



^1 



(71) 



This mass matrix is diagonalized by a unitary matrix T, 



TMfT-^ 



ml 








ml 



(72) 



In this scenario, the quantities ( ^ij) take a particularly simple form: 



VjlTal — \Vj2Ta2 

\uj2Tal 



for I 



b, s; a 



''1,2 



(73) 



while the sum over the squarks of the first two generations is determined by 
31. 



Let us for the moment neglect the mixing between and i^, i.e. consider 
the case mf^^ = 0, T = 1. Evaluating the first line of (|69D to lowest order, 
we find for the contribution of a light chargino and after summing over the 
different squarks: 



{mi 



1 



mj m 



w 



2 mb 



m?, / 



+0{al) + (subl.) +0(cot2/?) 



with the abbreviation 



~ 2 
mt 



~ 2 



m 



+ 



(74) 
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Similar, although rather lengthy expressions are obtained if the mixing be- 
tween and is taken into account, and analogous results are found for 
the other coefficients in (^) and (|70|) . As has already been pointed out in 



2^, the sign of the product Uj2Vji depends on the sign of /Xh, so that this 
leading contribution for large tan (3 can have either sign. 

A closer look at (0) shows two counteracting effects: a reduction of the 
leading coefficient due to QCD running from rhu down to rhj, while the 
term in square brackets shows an enhancement due to a "Q CD-softening" 
of the GIM cancellation, independent on whether is larger or smaller 
than rhu- The actual size of the corrections depends of course on the mass 
splitting between the squarks as well as on the splitting between the mass of 
the chargino and the squarks; since squarks can be an order or magnitude 
heavier than the lightest chargino, we estimate this coefficient to be of the 
order of 

14 a^i^rhj 

3 477 \^ rhj ' rht 

Similar results are found when analyzing the other expressions, so we will 
in general expect corrections up to 0(15%) with either sign. An exceptional 
situation occurs when, due to these super-GIM cancellations, the lowest- 
order contribution to O'l^ is accidentally lower than to Of^^ by orders of 
magnitude, since the above reasoning did not take into account the mixing 
of Olr into for scales below the heavy thresholds. In this sensible 
answer is obtained only when using the full expressions. 



In^, In^ <15% 



6 Conclusions 

We have extended the calculation of the leading QCD corrections for the 
inclusive b s'j decay to the MSSM in the framework of effective field 
theories. It was shown that it is important to properly treat the high-energy 
scale at which the particles in the loop are integrated out, as well as how to 
calculate the QCD corrections between if the masses of the particles in the 
loop are vastly different. To this end, we have calculated the leading order 
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anomalous dimension matrices for the operators for the various scenarios that 
are relevant to this process in the MSSM. 

We found that, while the SM contribution to the Wilson coefficients at the 
weak scale gets enhanced in the limit of a heavy top quark by about 15-20%, 
the contribution from a loop with a charged Higgs gets actually slightly 
reduced by a few percent. For the contribution from the chargino loops 
the result depends strongly on the mass spectrum of the squarks and the 
charginos as well as on the mixing angles. Typically, one expects corrections 
up to the order of 15% with either sign, which is however less than the 
enhancement of the SM contribution. 

Given a range of values for the inclusive decay, if one applies the above 
results to a parameter space analysis for a particular SUSY model, one will 
essentially find a relaxation of the bounds on the mass of the charged Higgs, 
especially in the region of large tan (3. The impact of the modification of the 
QCD corrections for the chargino loop contribution is not seen so easily, but 
we expect a smooth deformation of contours in analyses like [0, with 
the strongest effect in those regions where the lowest order contribution to 
the coefficient C'ogj^ is small although the chargino is relatively light. 

Finally we would like to point out that for the inclusive decay rate, even after 
taking into account the real gluon emission and virtual corrections below the 
h scale (03 , the leading order prediction remains uncertain by about 25% due 



to the residual scale dependence alone |]T0|, (for the amplitude it is of the 
order 10-15%). Once a full next-to- leading order calculation is available for 
the SM, it may be combined with the above results to obtain predictions in 
the MSSM with comparable precision. 
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A Wilson coefficients at one loop 

We quote here the results for the Wilson coefficients at one-loop order when 
both particles in the loop are integrated out at a common scale. These results 
will be used for the determination of subleading terms. They also provide 
an important cross-check for the leading terms obtained by the calculation 
in the effective theory, as well as for some of the entries in the anomalous 
dimension matrix. 

We find it convenient to use the following functions that appear in the eval- 
uation of the coefficients of the basis operators: 



Fi{x) 
F2(x) 



X — 5x — 2 xlnx 



12(x-l)3 2{x-iy 
+ 5x — 1 x'^lnx 
12(x - 1)3 ~ 2(x- 1)' 



^ , , X — 3 In X 



2(2; -1)2 (a;-l)3 
^ , , x + 1 X In X 



These functions are identical with those given in the appendix of ref. [|17 
Some of their properties are: 

Fi (-\ = xF2(x), F2 (-) = xFi(x), F4 (-) = xF4(x) 
Vx/ Vx/ Vx/ 

Fi(x) + F2(x) = ^F4(x) = ^ - ^xF3(x) 
xFi(x) + F2(x) = ^ 

F3 (-) = -X (F3(x) + 2F4(x)) + ^ (76) 
Vx/ X — 1 
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A.l Standard Model loop contributions 

Integrating out the W, the charged would-be Goldstone bosons and an up- 
type quark simultaneously, we obtain the one-loop expression of the Wilson 
coefficients of the effective Hamiltonian (|l|): 

Cqi = —xFiix) 



LR 



Cpl.i = Cpi.3 = hx + 2) (2F2(x) + Fs{x) + 2F^{x)) 

L L ,3 

Cpi.2 = hx + 2){F2{x)-F^{x)~2F^{x)) 

CplA = (x-2)F4(x) (77) 

Cp4 = (3 - 2x^3 (x) - 5x^4(0;)) 
^ Qh 

Here x = {mq/m^j^Y . Note that for large x all coefficient functions are 
bounded, except for Cp2, which grows logarithmically with x. For small x, 
Cpi.i, Cpi.2, Cpi.3 and Cp2 diverge logarithmically. 



A. 2 Charged Higgs loop contributions 

Integrating out the charged Higgs and an up- type quark simultaneously, the 
corresponding expressions are {y = {mq/mui)^): 

Col, = y^^(y^ 

CplA = CpUs = l-y{2F2{y) + Fs{y) + 2F,{y))cot'(3 

L ^ O 

Cpia = lyiF2iy)-F,iy)-2F,{y))coep 

CplA = -Q^Cp, = yF^{y)coe [5 (78) 
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cot^^ 



A. 3 Chargino loop contributions 



Finally we give the expressions for integrating out a chargino and an up- 



type squark. Setting z 



(m 



■,/maY-, where mj and represent the mass 



of the chargino Xj ^"^^ of the up-type squark Ua respectively, and using the 
couplings defined in eq. (pTI), one finds 



03 



C pi,i 



c 
c 



Pi 



mi, \ma J 


mb Krha J Qi 



-2)F^{z) 



Cpl,3 



g*jasgjah 



nin 













^F2{z) 










z 





g*jasgjah . "^W ^ 





\maJ 













2F2{z) - 










Z 





(79) 



g*jasgjah _ "^W 



After application of the equations of motion, these expressions are consistent 



with the corresponding expressions in |T7 
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